Chiral-mediated entanglement in an Aharonov-Bohm ring 
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We study the orbital entanglement in a biased Aharonov-Bohm ring connected in a four-terminal 
setup. We find that the concurrence achieves a maximum when the magnetic flux "J>s coincides with 
an integer number of half a flux quantum $o/2. We show that this behavior is a consequence of the 
existence of degenerate states of the ring having opposite chirality. We also analyze the behavior of 
the noise as a function of "1? and discuss the reliability of this quantity as an entanglement witness. 



(CME) because its creation is possible by the existence 
of intermediate states in the AB-ring which are coher- 
ent superpositions of two different chiralities for the elec- 
tronic motion. Remarkably, this kind of entanglement 
can be also defined in a transport setup, which is identi- 
cal to the one that was proposed to define the electron- 
hole entanglement 1 1 1 2 1 1 3 . In fact, we also show here 
that the noise current-current correlation function can be 
a good witness for the entanglement in our setup. 

The paper is organized as follows. In Sec. II, we 
present the model in detail and we sketch the formulation 
of the theory, where we calculate the reduced density- 
matrix of a two-particle system. The results of this work 
are presented in Sec. III. The conclusions are presented 
in Sec. IV. 



I. INTRODUCTION 

The increasing interest in quantum information pro- 
cessing is boosting the search for mechanisms to produce 
and control entanglement in devices of different nature^ 
Photonic devices are routinely used for preparing and de- 
tecting entangled photon pairs.— The main limitations of 
such devices is the non-existence of deterministic sources 
as well as the difficulty in controlling the interaction be- 
tween the photons. Entanglement mechanisms have also 
been proposed in solid state devices, like quantum dots, 3 
and superconductors p3 The fundamental ingredient be- 
hind all these mechanisms is a many-body interaction. 
More recently, it was determined that entanglement is 
also possible in systems of non-interacting electrons.^ 
In particular, it was shown that electron-hole entangled 
pairs can be produced by biasing a tunneling barrier. 

The edge states of systems in the quantum Hall regime 
can be employed in solid state devices to produce elec- 
tron beams with properties similar to those of a photon 
beam in optical setups. Exploiting this analogy, several 
theoretical and experimental proposals of electronic in- 
terferometers have been reported^ Interestingly, the elec- 
tronic counterpart of the Hanbury-Brown-Twiss device 
has been analyzed in a configuration of edge Hall states 
that do not have interfering orbits4^ For this reason, 
Aharonov-Bohm (AB) effect takes place only at the level 
of two-particle correlation functions, while the single- 
particle AB effect is not present. This entanglement 
seems to be related with an asymmetry in the device 
which favors the internal production of particle-hole pairs 
and manifest itself in the behavior of the current-current 
correlation functions. The relevance of the current- 
current correlation functions as entanglement witness has 
been also discussed in other fermionic systems 

The aim of this work is to establish the existence of or- 
bital entanglement in AB systems where single-particle 
interference is present. Using a microscopic model, we 
show that when coupling the AB-ring to four leads (two 
on the right and two on the left), the post select two- 
particle states of electrons at opposite leads are typically 
entangled. We name it "chiral-mediated entanglement" 



II. THEORETICAL TREATMENT 
A. Model 

We consider the AB single-channel ring with a 
magnetic flux. The ring is connected to four 
terminals^^ ' 10 ' 12 as shown in Fig. [TJ Two of the ter- 
minals, those labeled by a = 1, 2, at the left side of the 
ring, are at a higher voltage V with respect to the ones 
at the right, labeled by a = 3,4. All terminals are or- 
dinary single-channel metallic leads where electrons can 
move either to the right or to the left. For simplicity, we 
consider spinless electrons and we describe the setup by 
the Hamiltonian 



a=l 



H c ,a) + H. 



ring > 



(1) 



where H a = J^ fe £k a c\ Ck a are Hamiltonians of non- 
interacting electrons representing the leads. For the AB- 
ring we use a non-interacting model, where electrons 
move with velocity v either clockwise (+ chirality) or 
anticlockwise (— chirality). The Hamiltonian is 



H rmg = ^2 dxv\^l(x)V x ^ x {x), 



(2) 



being A the chirality, V x = — id x — 4>, 4> = $/(i$o), 
with <E> = 27T$s, where $b is the magnetic flux, $o = 



2 



hc/e is the flux quantum and L is the length of the 
ring. The contacts between the leads and the ring 
are modeled by tunneling terms of the form H ca = 
12k A=± Wk <* [ c fc 1> ^^( x a) + H.c.}, where x a define the po- 
sitions of the ring to which the leads are attached. We 
consider the leads to be at zero temperature. The chem- 
ical potentials enforce a bias voltage between left and 
right leads, i.e. p,\ = p 2 = Pl] pz = Pi = Pr and 
fJ>L—fi>R = eV, which we assume to be very small eV ~ 0. 



right or in the left leads ; 6 ' 10 ' 12 We introduce the operators 

-fin 



it — J J /it — J J 4! — J J /it = J J 



*00 — %i °fc 3 ' ^01 — ^kik,4> -"-10 — "k 2 Afc' H ~~ k 2 ^k 4 ' 

which create one particle in one of the left leads and a 
second particle at one of the right leads. The ensuing 
the 4x4 density matrix describes a system of two qubits 
with elements 

[ P {2) {e)U. a , w =J Lm25(e-e ka ){Ai b A a , b ,), 



B. Reduced Density Matrix and Concurrence 

In a setup as the one in Fig. [1] electrons tunnel from 
the left leads to the right ones. We aim to define the ef- 
fective density matrix describing the quantum state post- 
selected from the total two-electron state by projecting 
out the components where both electrons are either in the 



where Ao is a normalization factor, while the mean value 
is taken in the nonequilibrium state with a net current 
Sowing from the left to the right. We remark that in our 
approach matrix elements of are obtained in terms of 
operators appearing in the Hamiltonian H . Expectation 
values of four time dependent creation and/or annihi- 
lation operators are computed using the nonequilibrium 
Green function formalism and Wick theorem, 15 



(4jt)cUt)c kx (t)c ks (t)) = 



de 



de' 



G 



k x ,k a ( £ £ ) < ~"k !i ,k l3 



(3) 



where k (e) is the Fourier transform with respect to 



k a ,k 

t — t'oi the lesser Green function 

Gt, k At-t')=i(cUt')c ka (t)). 



(4) 



We assume that eV ~ and we are interested in analyz- 
ing e r~j fx = (p L + pr)/2. We present the corresponding 
expression of p^ in Appendix A. 

From the above density matrix we compute the con- 
currence, which is a good measure of entanglement^ In 
this case it is given by 



C 



(2) 



max 



{»• 



A, - 



where A.; are the eigenvalues of 



R 



(5) 



It 



P (7 y^ a vP a y® a y m decreasing order, 
interesting to notice that the concurrence so calculated 
is equivalent to the one obtained by the spin-dependent 
scattering formalism for a single mode conductor j££ 
which is 

C = 2 ^ ri(1 ' ri)T2(1 - T2) . (6) 

T\ + T 2 - TlT 2 

Here, t\,ti are transmission eigenvalues of the scatter- 
ing matrices s a ^{p) = 6 a< p - iy/TjT^G^(p), where 
Gj a (p) is the Fourier transform with respect to t — t' of 
the retarded Green function 



t') = -i@(t 



0^({*A(x a ,i),4(^/)}) 
A, A' 

(7) 

evaluated at p, being x a the position of the ring at which 



the wire a is attached, while r a = 27r^ A 

Sk a )>& 
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C. Noise 

As pointed out in Refs. HI andli~3l the concurrence can 
be expressed in terms of correlators that quantify the 
degree of violation of Bell inequalities. In transport se- 
tups the latter can in turn be directly related to current- 
current correlation functions, which are amenable to be 
experimentally detected. We, thus, turn to analyze the 
connection between these correlation functions and the 
above discussed entanglement. We first outline the pro- 
cedure to compute the current-current noise within our 
treatment. The current passing through the contact to 
the terminal a, can be expressed by the following oper- 
ator J a (t) = (e/h)Y, ka w ka [icl a (t)^(x a ,t) + h.c]. The 
zero frequency shot-noise is a measure of the current- 
current correlations in different terminals. It reads 

<W(0) = i J dr({SJ a (T),SJ p (0)}), (8) 

where 5J a = J a — ( J a ). We calculate (ec.©) by evaluat- 
ing a bubble diagram in terms of nonequilibrium Green 
functions. The corresponding expressions are presented 
in Appendix IIVI 



III. RESULTS 

A. Qualitative Analysis 

To understand the origin of CME it is useful to be- 
gin analyzing the electronic states of the isolated AB- 
ring. The Hamiltonian can be diagonalized in momentum 



3 



space: Denning ^x(x) — l~/VN^2 p e %px c Pt x with Af a 
normalization factor, p — 2-KnjL and n £ Z we obtain 

H rmg = J2x Ep £ p,A(*)Cp,A c p,>' with £ p+ = *v(p-4>) ( a 

cutoff in the single- particle energy spectrum is assumed). 
The effect of the magnetic flux on the energies e P) a($) is 
illustrated in Fig. [T] Depending on the magnetic flux, 
there can be zero, one or two single-particle states \p, A) 
with a given energy. States with different chirality A = ± 
are degenerate only when the flux is an integer multiple 

Of 7T$ . 

Let us consider first the case where two degenerate 
states with opposite chiralities exist in the ring. These 
two states behaves as an intermediate qubit that couples 
to the qubits defined by the leads. We will argue that 
in this case CME naturally emerges between electrons 
at the right and left leads. The iV-particle states with a 
Fermi energy e p can be obtained from |0) , that represents 
the Fermi sea with N — 2 particles filling the states with 
£ p,a(^) < e F , where v(pp — (/))= ep. Thus, we have 
luring) = c PF i+ c PFi _ |0) . When the ring is in contact 
with the four leads, particles can tunnel between the ring 
and the reservoirs. For weak coupling we can assume that 
each of the chiral levels with p p hybridize with the levels 
of the leads having the same energy e p . This is described 
by the following effective Hamiltonian 

He i} =wJ2 E e lXp ^[clc PF , x + H.c.}, (9) 

a=l A=± 

where w is the effective tunneling parameter, c\ creates 
an electron in the single-particle state of the a-lead with 
energy £f. a = £f (we take ep = without loss of gen- 
erality). This Hamiltonian has four eigenstates of the 
form 

4 

\ipn) = E 7n ' aC L + E 7 ™' aC L,a]I°>' ra = l,..., 4, 

a = l A=± 

(10) 

where the coefficients 7„ are the weights of the eigen- 
states in the chosen base. It also has two additional de- 
generate states of the form \ip n ) — Yl a =i "K" c ! n = 
5,6. The latter correspond to states that do not hy- 
bridize with the ring. When two particles are present, 
two different such states must be occupied. It is simple 
to show that any state of this type has a sizable projec- 
tion on states of the form J2 a ^/3 ^■<x,P c k a c k0 1^) f° r some 
non- vanishing coefficients A Qj/ 3. This two-particle state is 
typically entangled in the orbital indices a, ft of opposite 
leads. Notice that it is also possible to use two AB-rings 
with degenerate levels as a two-qubit system coupled by 
two conducting leads (intermediate qubit). 

On the other hand if we consider the case of a one 
chirality ring (e p = v(p — 4>) for example) the situation 
drastically changes and no significant entanglement be- 
tween left and right leads is attained. This is can be 
seen because the effective Hamiltonian has a different 
level structure. In fact, H e ff can be naturally written 




(b) 

i cp/<t> =a% <p/q> *n?c 

V V 

FIG. 1. (a) Scheme of an Aharonov-Bohm ring with two 
chiralities, attached to four leads. The chemical potentials 
are such that there is a bias voltage between the left and right 
leads. Electrons incoming from the left can tunnel to the ring 
and escape through leads, (b) Linear dispersion relation of an 
isolated AB-ring, where electrons move with velocity v either 
clockwise (+ chirality) or anticlockwise (— chirality). The 
magnetic flux determines whenever the Fermi level is two- 
degenerate or not. 

in terms of operators that are linear combinations of the 
lead operators c' k as H e ff = w{f^c PFt \ + H.c), with 

P = (1/2) J2a=i e %PFXa c { k , while there are three addi- 
tional orthogonal linear combinations of these operators 
which do not hybridize with the ring. The two eigen- 
states of H e f f are linear combinations of a single-particle 
state of the leads and a single-particle state of the ring. 
Thus, a two-particle state of this Hamiltonian has never 
the two particles in the leads and therefore no orbital 
entanglement is possible. 

The above argument suggests that by varying the mag- 
netic flux, or the chemical potential of the leads (or equiv- 
alcntly, a gate voltage applied at the ring) we can induce 
the system to switch from a situation with no orbital 
entanglement between the leads onto another situation 
with inter-lead entanglement. This is done by varying 
$ and/or /i in order to have degenerate chiral states of 
the ring at the Fermi energy. We now present a rigor- 
ous calculation of the entanglement for the states that 
are relevant for a transport experiment in the coherent 
regime and explain how the CME depends on flux and 
chemical potential. We also discuss the way in which it 
may be detected in transport experiments. 

B. Numerical Results 

1. Concurrence 

The behavior of C as a function of the mean chemical 
potential /_t and the magnetic flux $ is shown in Fig. [2] 
The concurrence is maximal at $/$o — n7T and for values 
H close to the energy of two degenerate chiral states of 
the ring. The same type of behavior is observed for other 
configurations of the wires, corresponding to contacts at 
different positions x a . For this Fig. we considered wires 
with a bandwidth W a and T a = il a y/W% — li 2 , where 
Q a is a constant, but the same behavior is obtained for 
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Jt/2 n 3k/2 2% 



<3>/<E> 

FIG. 2. (Color online) Concurrence as a function of &/&o 
and mean chemical potential /i for a ring with length L = 20 
with wires connected at xi = 1, x% = 6, £3 = 11, £4 = 16. 




7i/2 ji 37i/2 2ji ji/2 ji 37t/2 2% 



<!>M> '^/'I'q 

FIG. 4. (Color online) Shot noise (squares) between left and 
right terminals, as a function of $/$o for two values of the 
chemical potentials at which there are degenerate states, (a) 
/i = 0.3175, (b) /1 = 0.7866. The corresponding plots for the 
concurrence are also shown for comparison (circles). 



2. Noise 
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FIG. 3. (Color online) Concurrence as a function of $/$o 
and mean chemical potential /i for a ring with only one chi- 
rality. The ring parameters are the same as in Fig[2] 



leads with a constant density of states. For some chemical 
potentials C exhibits maxima at $/$o = 0, mod (2ir), 
which corresponds to the energy of two degenerate states 
of the ring, but achieves again the maximum value within 
a wide range of fluxes centered at <!>/<I>o = 7r, mod (2tt). 
This feature is analyzed below in more detail. 

Instead, if we evaluate C for the Hamiltonian @ re- 
stricted to a single chirality, we find negligibly orbital 
entanglement in the leads within the whole range of <£> 
and li. These results are presented in Fig|3J In this case 
the ring behaves as a single-level system, and prevent the 
formation of orbital entangled states at the leads. 



We now turn to present results on the behavior of the 
current-current correlation functions. Our aim is to an- 
alyze if the signatures of entanglement found in the be- 
havior of the concurrence can be also identified in the be- 
havior of the noise. Results for the total left-right noise 
correlations S = Y] a —i 2 s=3 4^a,/3(0) (equal to the self- 
correlation — J2 a o =1 2 S a ^(0) ) are shown in Fig. |4j 
The left (right) panel, corresponds to a chemical poten- 
tial \i for which there are two degenerate chiral states in 
the ring for <&/^>o = 0, (71"), mod(27r). The behavior of 
the concurrence is also plotted for comparison. In both 
cases, S, along with C exhibit maxima at the fluxes for 
which the two degenerate chiral states are resonant at 
the given ll. This points to the idea that noise is in- 
deed a reliable witness of orbital entanglement, as dis- 
cussed in the context of other electronic setups££&ii In 
the case of the left panel, both quantities are maximum 
at $/$o = 0,mod(27r). Within a range of fluxes which 
scans the width of the resonant degenerate levels of the 
ring they first decrease and then increase, displaying a 
dip. As the flux increases further, the behavior of these 
two quantities, however, depart one another. While S 
tends to vanish around < i > / < E > o = i") the concurrence dis- 
plays a wide plateau with height C ~ 1. Qualitatively, 
the same type of behavior is observed in the right panel. 
In this case, both quantities exhibit a sharp maxima at 
resonance (see the peaks around <!>/<f>o = 7r). The concur- 
rence displays a plateau and another (lower) maximum 
around ^/^o = 0,mod(27r) while S is vanishing small. 
On general grounds this is rather surprising, since one 
could easily imagine situations with a sizable noise with- 
out entanglement, but here we have the converse situa- 
tion. 

In order to further understand the relation between 
the behavior of C and S as well as the connection to the 
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FIG. 5. (Color online) Probability of finding two particles at 
the left and right leads as a function of magnetic flux <I> / <E>o . 
Panels (a) and (b) correspond to the same parameters of Fig. 
[4] (a) and (b), respectively . 



CME, we analyze the probabilities 

a k a 

p R r = UJ2 S ( £ - e ^)( c lc k3 ci^) (ii) 

of finding two particles at the left (right) leads, respec- 
tively. In Fig. [3] we show the behavior of these two 
quantities for the same parameters of Fig. 2J It is 
clear that Pll and Prr are both sizable when tunneling 
through the ring is allowed. For these chemical poten- 
tials, this corresponds to the narrow window of fluxes 
around < J>/'I > o = 0(7r), for the case of the left (right) 
panels, respectively, within which the degenerate levels 
of the ring remain resonant. Beyond these values, the 
transmission from left to right is blocked and the two 
particles have a high probability of remaining within the 
left wires. A further analysis of Figs. [5] and 0] in the 
light of the results shown in Fig. [5j then reveals that 
the high values of concurrence in the plateaus away from 
$/$o = 0(7r), mod(27r) in the case of the left (right) pan- 
els of Fig. SI correspond to states which have a very low 
probability of taking place. Instead, the resonant situa- 
tion with two degenerate chiral states of the ring at the 
Fermi energy, leads to a high orbital entanglement which 
is clearly witnessed by a high noise amplitude S. 



IV. CONCLUSION 

To summarize, we introduced a new mechanism for 
orbital entanglement. This type of entanglement is origi- 
nated in the spectral nature of the AB ring and is highly 
sensitive to the magnetic field. Orbitally entangled elec- 
tronic pairs can be produced by suitably tuning the mag- 
netic field, the chemical potential or a voltage gate at 
the ring, in order to have a degenerate pair of electronic 
states with different chiralities at the Fermi energy (inter- 
mediate qubit). In fact, if the right and left leads are me- 
diated by a single-level system (single-chirality ring) the 
orbital entanglement is negligible. This type of entangle- 
ment can be detected in transport experiments, the shot 
noise being a good witness. The setup of Fig. Q] could be 
experimentally realized in an architecture based on the 
quantum Hall regime of a 2D electron gas, by substitut- 
ing each wire by a pair of incoming and outgoing edge 
states and the ring by a pair of edge states with different 
chiralities, separated by a narrow circular wall. In such 
a setup it would be possible to combine the main system 
of Fig. [T]with beam splitters connected at the wires, in 
order to test Bell inequalities and even perform the full 
quantum tomography following the protocol of Ref. ITol . 
Two interesting possible generalizations are the combina- 
tion between static and flying qubits by combining with 
quantum dotsj^i as well as the introduction of dynamical 
single-particle emitters at the sources^ 
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APPENDIX 

A. Reduced Density Matrix 

We evaluate the post-selected state of two electrons 
at opposite leads, p^ 2 \ in terms of Green functions (the 
Fourier transform of ec.©). The explicit expression of a 
general matrix element, up to a normalization factor, is: 
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II E<^ - *«M a (*)<*»(*)<%,(*)} = T 2 [/ a ( £ )/^( £ )Gf iQ ( £ )G^( £ ) - Ue)f s (e)G?Js)Gfjs) 



4 

- h (e)Me)G?Js)Gfje)]- l J2r 3 [Me)Ms)G?Je)G?Je)G*Je) 



7=1 



f a (e)f y (e)G*Je)G?Je)G$Jej\ + £ ^f v (s)Me)Gfje)Gfje)G^(e)G^Je) 



7 ,7J = 1 



+r 



i/ a ( e )/ fl ( e )G» / ,( e )+t/ a (e)/,( e )G^(e)-r£;/ 7 ( e )/ (e)G^( e )G^( e ) 



7=1 
4 



^( £ )/ Q ( £ )Gf :Q ( £ ) + l / A ( £ )/ 5 ( £ )G^(e)-r^/ 7 ( e )/ 5 ( £ )Gf !7 ( e )G^( £ ) 



7=1 



h s ,kf) 



hot &5 



(12) 



where we considered all spectral densities equal T a = T Note that G^Je) — J2\ \> Gf" x , (x a , Xp, e).The retarded 
Vcc. / a (e) is de Fermi distribution function of the lead a. Green functions are evaluated by solving the Dyson equa- 



tion, 



G£ v (x, a/, e) = <?f, (x, x', e)S(x - x')5 x ,x> + E E G £a» (x, b 7 , e)E?(e)$£ (x 7 , a/, e) 

7=1 A" 



(13) 



being 



of reservoir a. 



fea 



sj(a; ) i / ,e) = — ^ 



-zfc(a;— a;') 



X , ^ e- £fc.A'( $ ) + V 

fc — — Kn 



(14) 



B. Shot Noise Calculation 



where fco is the energy cut off and A4 = 2kg + 1 a normal- 
ization factor, the uncoupled retarded Green function of Here we show the shot noise expression of ec.([S]) in 
the ring and £ 7 (e) = Ylk \ w k a \ 2 9k ( e ) the self energy terms of retarded Green functions. 

I 



15 



(15) 



e 

2ft 



^ L E [r«I>r 4 r 7 [(l - m)fs{e) + (1 - .f A -( £ ))/ 7 ( £ )]G« 7 ( £ )G* ( £ )Gf,,( £ )G* ( e )] 



"5,7=1 



-2r a I>[(l - fp(e))f a (e) + (1 - /a(e))/fl(«G^(e)G?,«(E)} 

4 

-2r a r^e{iGf )a ( £ ) X] Gf; 5 (e)G^( £ )r 5 [(1 - f a (e))f s (s) + (1 - /*(e))/a(e)]} 

5=1 
4 

-2r a r^e{iG«^( £ ) Gf, s (e)G% s (e)r s [(1 - /^(e))/s(e) + (1 - /ife))//,^)]} 

5=1 

4 

-2<5 Q ^r Q J2 r5|G^ 5 (e)| 2 [(1 - fs(e))f a (e) + (1 - f a (e))fs(e)} . 

5=1 ' 
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